‘CROCOPY  RESOLUTION  TEST  CHART 

national  bureau  of  STANDARDS-1963-A 


PTrr  rnpT  ad-ai7i  434 


AFWAL-TR-86-2004 
Volume  I 


COMPUTATIONAL  FLUID  DYNAMIC  STUDIES  OF  CERTAIN  DUCTED 
BLUFF-BODY  FLOWFIELDS  RELEVANT  TO  TURBOJET  COMBUSTORS 

VOLUME  I  -  Time-Dependent  Calculations  with  the  k-e 
Turbulence  Model  for  an  Existing  Centerbody  Combustor 


M.  S.  Raju  and  L.  Krlshnamurthy 
University  of  Dayton 
Research  Institute 
Dayton,  Ohio  45469-0001 


July  1986 


FINAL  REPORT  FOR  PERIOD  16  AUGUST  1984  -  30  SEPTEMBER  1985 


Approved  for  Public  Release,  Distribution  Unlimited 


AERO  PROPULSION  LABORATORY 

AIR  FORCE  WRIGHT  AERONAUTICAL  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 

HRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-6563 


DTIC 


SELECTE 
AUG  2  7 1986 


B 


HOTICE 


"When  Government  drawings,  specifications,  or  other  data  are 
used  for  any  purpose  other  than  in  connection  with  a  definitely 
related  government  procurement  operation,  the  United  States 
Government  thereby  Incurs  no  responsibility  nor  any  obligation 
whatsoever;  and  the  fact  that  the  Government  may  have  formulated, 
furnished,  or  In  any  other  way  supplied  the  said  drawings, 
specifications,  or  other  data.  Is  not  to  be  regarded  by 
Implication  or  otherwise  as  In  any  manner  licensing  the  holder  or 
any  other  person  or  corporation,  or  conveying  any  rights  or 
permission  to  manufacture,  use  or  sell  any  patented  invention 
that  may  in  any  way  be  related  thereto.” 

This  report  has  been  reviewed  by  the  Office  of  Public  Affairs 
(ASD/PA)  and  Is  releasable  to  the  National  Technical  Information 
Service  (NTZS) .  At  HTIS,  It  will  be  available  to  the  general 
public.  Including  foreign  nations. 

This  technical  report  has  been  reviewed  and  Is  approved  for 
publication. 


UJ.  rn.  _ 

W.M.  ROQUEMoS 
Fuels  Branch 

Fuels  and  Lubrication  Division 
Aero  Propulsion  Laboratory 


ARTHUR  V.  CHURCHILL. 

Chief,  Fuels  Branch 

Fuels  and  Lubrication  Division 

Aero  Propulsion  Laboratory 


FOR  THE  COMMANDER 


ROBERT  D.  SHERfRILL,  Chief 


Fuels  and  Lubrication  Division 
Aero  Propulsion  Laboratory 


"If  your  address  has  changed.  If  you  wish  to  be  removed  from 
our  mailing  list,  or  If  the  addressee  Is  no  longer  employed  by 
your  organization,  please  notify  AFWAL/POSF,  W-PAFB,  OH 
45433-6563  to  help  us  maintain  a  current  mailing  list. 


Copies  of  this  report  should  not  be  returned  unless  return  Is 
required  by  security  considerations,  contractual  obligations,  or 
notice  of  a  specific  document. 


Unclassified 


•CCUAITV  CLASSIFICATION  OF  THIS  PAGE 


la.  REPORT  security  CLASSIFICATION 


/in<  W3V 


REPORT  DOCUMENTATION  PAGE 


2».  SECURITY  CLASSIFICATION  AUTHORITY 

N/A 


2b.  OECLASSIFICATION/OOWNGRAOING  SCHEDULE 

N/A 


4.  PERFORMING  ORGANIZATION  REPORT  NUMBERISI 

UOR- rR-85-82 


Sa.  NAME  OF  PERFORMING  ORGANIZATION 

University  of  Dayton 
Research  Institute 


6€.  AOORESS  (City,  State  and  ZIP  Cod€) 


Dayton,  OH  45469-0001 


.  NAME  OF  FUNOING/SPONSORING 
ORGANIZATION 


8e.  AOORESS  (City.  Statt  and  ZIP  Coda) 


1b.  RESTRICTIVE  MARKINGS 

N/A 


3.  OISTRIBUTION/AVAILABILITY  OF  REPORT 

Approved  for  Public  Release,  Distribution 
Unlimited 


5.  MONITORING  ORGANIZATION  REPORT  NUM6ER(S) 

AFWAL-TR-86-2004 .  Volume  I 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Aero  Propulsion  Laboratory  (AFSC) 
Air  Force  Wrioht  Aeronaut! 


7b.  AOORESS  Wity.  Slate  and  ZIP  Code) 

AFWAL/POSF 

Wright-Patterson  Air  Force  Base,  OH 


8b.  OFFICE  SYMBOL  9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

Ilf  applicable,  F33615-84-C-2411 ,  Task  02 


11.  title  (fncludt  Security  Clamificationt 

id 


12.  PERSONAL  AUTHORISI 

M.  S.  Raju  and  L.  Krishnamurth 


13a  TYPE  OF  REPORT 

Final  Report 


1 10.  SOURCE  OF  FUNDING  NOS. 

PROGRAM 

PROJECT 

TASK 

WORK  UNIT 

ELEMENT  NO. 

NO. 

NO. 

NO. 

62203F 

3048 

05 

43 

13b.  TIME  COVERED 
FROM 


14.  DATE  OF  REPORT  (Yr,  Mo.,  Day) 


TO  9/-30/85  July  1986 


1B.  SUPPLEMENTARY  notation 

Work  performed  under  the  Scholarly  Research  Program 


16.  PAGE  COUNT 

46 


17. 


FIELD 


COSATI  COOES 


IB.  SUBJECT  TERMS  iContinue  on  rwvtrte  if  necrcaory  and  identify  by  block  number) 

Bluff-Body  Near  Wake  Recirculating  Flows 

Confined  Turbulent  Flows  (See  Reverse  5idel 

Nnnreaf.t.inn  FIolic _  Keversc  iioe; 


19.  ABSTRACT  (Continue  on  reuerae  if  neceuary  and  identify  by  block  number) 

■'■>'4  numerical  investigation  of  the  near-wake  region  in  a  ducted  bluff-body  combustor 
by  finite-difference  computations  is  reported.  The  numerical  predictions  are  based  upon 
the  time-dependent,  compressible  Navier-Stokes  equations  and  the  k-e  turbulence  model. 

The  standard  k-_e  turbulence  model  -has-beert^^ihBdified  to  account  for  the  nonstationary 
terms.  The  time-dependent  calculations ‘Ttev^ addressed  the  nonreacting  near-wake  flowfield 
of  the  centerbody  combustor  with  only  the  annular  air  stream  present. 

Flowfield  predictions  for  a  com^stor  inlet  mass  flow  of  2  kg/s  with  the  time-depend¬ 
ent  formulation  incorporating  the  kniC  turbulence  model  show  the  attainment  of  a  steady- 
state  recirculating  flow  in  the  near  wake.  The  slow  axial  migration  of  the  recirculation 
vortex  towards  the  exit  boundary  which  was  noticed  in  the  earlier  time-dependent  calcu¬ 
lations  without  a  turbulence  model  is  no  longer  present.  Thus,  tli^present  results  have*^ 
eliminated  the  appearance  of  reverse  flow  \  (See  Reverse  Side)  *  1 

3a  OISTRIBUTION/AVAILABILITY  OF  ABSTRACT 
UNCLASSIFIEO/UNLIMITED  QO  SAME  AS  RPT.  □  OTIC  USERS  □ 


23a  name  of  responsible  individual 

W.  M.  Roquemore 


31.  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


aZb.  TELEPHONE  NUMBER 
(include  Area  Code) 

(513)  255-6813 


rrrTTTm 


33g.  OFFICE  SYMBOL 

AFWAL/POSF 


Unclassified _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


BLOCK  11  (Continued) 

Computational  Fluid  Dynamic  Studies  of  Certain  Ducted  Bluff-Body  Flowfields  Relevant  to 
Turbojet  Combustors  (U) 


Volume  I 

Time-Dependent  Calculations  with  the  k-e. Turbulence  Model  for  an  Existing  Centerbody 
Combustor 


BLOCK  18  (Continued) 

Unsteady  Flow 
Vortex  Shedding 


BLOCK  19  (Continued) 

at  the  exit  boundary  with  the  consequent  incompatibility  of  the  boundary  conditions,  and 
thereby  the  spurious  shedding-like  behavior  observed  previously.  The  steady-state  results 
in  the  present  study  dei-aonstrate  internal  consistency  with  the  time-averaged  measurements 
and  predictions  for  the  locations  of  the  vortex  center  and  the  centerline  rear  stagnation 
point.  Preliminary  computations  for  both  laminar-like  and  turbulent  flows  with  time- 
dependent  perturbations  of  inflow  boundary  conditions  do  not  show  any  oscillatory 
behavior  in  the  interior  of  the  combustor.  . 

Present  time-dependent  turbulent  computat+of^  with  the  MacCormack  algorithm  confirm 
the  earlier  observations  (from  the  laminar-like  calculations)  of  the  inability  of  this 
computational  procedure  to  accurately  capture  the  dynamic  features  of  the  ducted  center- 
body  combustor  flowfields.  It  is  likely  that  the  experimental  observations  of  the 
dynamic  behavior  of  this  configuration  in  reacting  flows  may  have  arisen  from  significant 
coupling  between  the  duct  acoustics  and  combustion  heat  releasei 


/ 


linclassified 


SECUniTV  CLASSIFICATION  OF  THIS  FAGS 


PREFACE 


This  final  report  was  submitted  by  the  University  of  Dayton, 
under  Contract  No.  F33615-84-C-2411 ,  Task  02.  The  research  task, 
entitled  "Assessment  of  Time-Dependent  Calculations  for  Gas 
Turbine  Combustor-Type  Flows,"  wais  sponsored  by  the  Scholarly 
Research  Program  of  the  Air  Force  Wright  Aeronautical 
Laboratories,  Aero  Propulsion  Laboratory,  Wright-Patterson  Air 
Force  Base,  Ohio,  under  Project  No.  3048,  Task  05,  Work  Unit  43. 
Dr.  William  M.  Roquemore,  AFWAL/POSF,  was  Project  Engineer.  The 
research  work  dealing  with  the  time-dependent  calculations 
reported  in  Volume  I  was  Initiated  in  August  1984  and  completed 
in  February  1985.  The  work  describing  the  time-averaged 
calculations  reported  in  Volume  II  was  initiated,  under  an 
Extension  of  Task  02,  in  April  1985  and  completed  in  August 
1985.  The  numerical  calculations  were  performed  by  Dr.  M.  S. 
Raju,  Associate  Research  Engineer,  and  the  report  was  written  by 
Dr.  L.  Krlshnamurthy ,  Senior  Research  Engineer.  The  authors 
express  their  appreciation  to  Dr.  Roquemore  for  his  interest  and 
to  Mr.  J.  S.  Stutrud,  AFWAL/POSF,  for  his  help  with  their  MODCOMP 
computer  in  the  time-averaged  calculations.  The  time-dependent 
computations  were  carried  out  on  the  CRAY-IS  computer  at  the 
NASA  Lewis  Research  Center  and  the  authors  express  their 
appreciation  to  Dr.  E.  j.  Mularz  and  Mr.  R.  W.  Claus  for  their 
support  of  this  computing  activity.  The  authors  also  acknowledge 
the  assistance  of  Ms.  Ellen  Bordewisch,  Ms.  Teresa  Harding,  and 
Ms.  Marlene  Haas,  UDRI,  in  preparing  this  report,  and  Ms.  Anita 
Cochran,  UDRI  in  technical  editing  thereof. 


DTIC 

Selecteb^ 
AUG2  7e88  J  I 

B 


ill 


TABLE  OF  CONTENTS 


SECTION  PAGE  NO. 

I .  INTRODUCTION  1 

1 .  Background  1 

2.  Implications  3 

3.  Previous  Work  4 

4.  Scope  of  Present  Work  7 

5.  Outline  of  Report  9 

II.  TIME-DEPENDENT  FORMULATION  10 

1.  Governing  Equations  10 

2.  Boundary  and  Initial  Conditions  13 

a.  Inflow  Boundary  Conditions  14 

b.  Outflow  Boundary  Conditions  14 

c.  Solid-Wall  Boundary  Conditions  14 

3.  Computational  Details  15 

a.  Lamlnar-Llke  Cases  with  Time-Dependent  Inflow  15 

(1)  Case  1  17 

(a)  Inflow  Conditions  17 

(b)  Outflow  Conditions  17 

(2)  Case  2  17 

(a)  Inflow  Conditions  17 

(b)  Outflow  Conditions  17 

(3)  Case  3  18 

(a)  Inflow  Conditions  18 

(b)  Outflow  Conditions  18 

b.  Turbulent  Case  with  Time-Dependent  Inflow  18 

III.  RESULTS  AND  DISCUSSION  19 

1.  Influence  of  Turbulence  Model  19 

2.  Effect  of  Time-Dependent  Inflow  Conditions  28 

a.  Laminar-Llke  Calculations  28 

b.  Turbulent  Flowfleld  Results  37 

IV.  CONCLUSIONS  AND  RECOMMENDATIONS  40 

1.  Conclusions  40 

2.  Recommendations  43 

REFERENCES  45 


V 


LIST  OF  ILLUSTRATIONS 


FIGURE 
1 . 

2(a)-(b) . 
2(c)-(d). 

3a. 

3b. 

4(a)-(b) . 

4(c)-(d) . 

5a. 

.‘ib. 

6a. 

6b. 


Baseline  Finite-Difference  Mesh  (60  x  46) . 


Velocity-Vector  Plots  for  Turbulent  Flow. 

Vortlclty-Contour  Plots  for  Turbulent  Plow. 

Velocity-Vector  Plots  for  Lamlnar-Lllce 
Calculations  (Case  1)  with  Forcing. 


Vorticity-Contour  Plots  for  Laminar-Like 
Calculations  (Case  1)  with  Forcing. 


Velocity-Vector  Plots  for  Leunlnar-Like 
Calculations  (Case  2)  with  Forcing. 


Vorticity-Contour  Plots  for  Laminar-Like 
Calculations  (Case  2)  with  Forcing. 


Velocity-Vector  Plots  for  Laminar-Like 
Calculations  (Case  3)  with  Forcing. 


Vorticity-Contour  Plots  for  Laminar-Like 
Calculations  (Case  3)  with  Forcing. 


Velocity-Vector  Plots  for  Turbulent  Plow 
with  Forcing. 


Vortlclty-Contour  Plots  for  Turbulent  Flow 
with  Forcing. 


PAGE  NO. 
16 

20-21 

22-23 

29 

30 

31-32 

33-34 

35 

36 

38 

39 


SECTION  I 
INTRODUCTION 


This  final  report  (in  two  volumes)  documents  the  results  of 
the  computational  fluid  dynaunlc  (CFD)  research  performed  by  the 
University  of  Dayton  for  the  Air  Force  Wright  Aeronautical 
Laboratories,  Aero  Propulsion  Laboratory  (AFWAL/PO) .  The  focus 
of  this  CFD  research  is  the  confined,  turbulent  recirculating 
flowfield  behind  a  bluff  body  in  the  AFWAL/POSF  research 
combustor . 

In  an  existing  version  of  the  POSF  combustor,  a  centerbody 
configuration.  Involving  the  turbulent  mixing  and  combustion  of 
an  annular  air  stream  and  a  central  fuel  Jet  in  the  near  wake 
region  of  an  axlsymmetrlc  bluff  body,  has  been  the  subject  of 
extensive  diagnostic  and  predictive  research.  Computations  of 
this  configuration  for  the  nonreacting  flow  due  to  the  annular 
air  alone  have  been  performed  with  the  time-dependent  Navler- 
Stokes  equations,  incorporating  a  realistic  turbulence  model, 
and  these  results  are  presented  here,  in  Volume  I  of  this  report. 
The  time-averaged  predictions  based  upon  a  solution  of  the 
Reynolds-Averaged  Navier-Stokes  equations  for  a  proposed  POSF 
combustor  configuration.  Involving  two  annular  air  streams — a 
swirl-free  outer  stream  and  a  swirling  inner  stream,  a  central 
fuel  Jet,  and  a  centerbody  imbedding  the  swirler  and  the  fuel 
nozzle,  are  reported  in  Volume  II. 

1 .  BACKGROUND 

The  performance  of  a  gas  turbine  is  significantly  Influenced 
by  the  fluid  mechanical  and  combustion  processes  in  the 
combustion  chamber.  The  details  of  the  fluid  motion  and 
associated  turbulence  largely  determine  the  mixing,  combustion, 
and  heat  transfer  characteristics  of  the  combustor.  Thus,  an 
accurate  prediction  of  turbulent  mixing  and  combustion  would 
require  the  knowledge  of  turbulence  structure  and  turbulence 
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scales  of  a  combustor  flowfleld.  The  viewpoint  adopted 
for  the  prediction  of  turbulent  flows  is  that  they  are  the 
unsteady  solutions  of  the  Navler-Stokes  equations^  and  that  their 
description  involves  probabilistic  features.  Such  a  viewpoint, 
however,  renders  the  prediction  of  turbulent  reacting  flows  in 
realistic  configurations  like  gas-turbine  combustors  exceedingly 
difficult,  since  the  full  unsteady  nature  of  the  flowfleld  must 
be  computed  even  if  one  is  interested  only  in  its  time-averaged 
behavior . 

In  principle,  one  can  obtain  the  time-averaged  information  in 
two  ways:  either  by  solving  the  time-dependent  Navier-Stokes 
equations  and  averaging  the  solutions,  or  by  averaging  these 
equations  first  and  then  solving  them.  In  practice,  however, 
both  avenues  suffer  from  major  difficulties.  Averaging  the 
Navier-Stokes  equations  results  in  equations  (whether  they  are 
the  Reynolds-averaged  equations  in  constant -density  flows,  or  the 
Favre-averaged  equations  in  variable-density  flows)  that  are 
confronted  by  the  Indeterminacy  known  as  the  closure  problem  and 
by  the  consequent  necessity  to  introduce  a  turbulence  model .  The 
direct  and  explicit  computation  of  the  time-dependent  equations, 
on  the  other  hand,  has  two  drawbacks.  First,  a  typical  turbulent 
flowfleld  is  characterized  by  a  large  number  of  Interacting 
length  scales,  whose  explicit  computation  by  the  numerical 
solution  of  the  full  equations  at  high  Reynolds  numbers  requires 
enormous  computing  effort.  The  second  source  of  difficulty  is 
the  apparent  lack  of  uniqueness  of  the  solutions:  if  any  two 
realizations  in  time,  with  apparently  identical  initial  and 
boundary  conditions,  are  distinctly  different,  averaging  makes 
sense  only  if  all  such  realizations  and  their  associated 
probabilities  are  known  (Reference  1). 


1  Although  the  fluid  mechanical  turbulence  is  generally 

accepted  to  be  contained  in  and  explained  by  the  Navier- 
Stokes  equations,  they  remain  to  be  fully  tested  for 
completeness  by  computations. 
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IMPLICATIONS 


2  . 


The  foregoing  difficulties  notwithstanding,  CFD  prediction  of 
complex  turbulent  flows  in  practical  systems  has  proceeded 
vigorously,  based  upon  computations  of  the  averaged  equations. 
Indeed,  such  computations,  involving  the  two-equation  model  (for 
the  turbulent  kinetic  energy  k  and  the  rate  of  its  dissipation 
e),  appear  to  provide  reasonable  predictions  of  the  overall 
trends  in  gas  turbine  combustor-type  flows  and  represent  the 
accepted  design  tool  in  the  industry  (Reference  2).  Neverthe¬ 
less,  we  must  recognize  that  these  time-averaged  computations 
are  based  upon  relatively  well-developed  statistical  approaches 
which  assume  the  classical  description  of  turbulence — a 
continuous  cascading  of  large  eddies  into  smaller  dissipative 
scales  with  random  interactions  thereamong.  Consequently,  the 
quantitative  predictions  have  acceptable  accuracy  only  when  we 
have  the  so-called  "fully  developed  turbulence."  Unfortunately, 
realistic  turbulent  flows  of  engineering  importance  do  not  become 
fully  developed,  and  the  predictive  success  therein  requires  the 
experimental  calibration  and  optimization  of  the  free  parameters 
(arising  from  the  closure  of  the  equations) . 

An  implication  of  greater  significance  for  predictors  is  the 
existence  of  large-scale  ("coherent")  structures,  which  has 
received  considerable  attention  recently  from  both  turbulence  and 
combustion  researchers.  Recent  experiments  (see,  e.g..  Reference 
3)  have  led  to  the  growing  realization  that  the  vortlcity 
fluctuations  in  turbulence  are  not  quite  so  random  or  incoherent 
as  was  previously  supposed.  This  viewpoint  suggests  that 
turbulence  should  be  computed  as  the  creation,  evolution, 
interaction,  and  decay  of  large-scale  vortex  structures. 

Although  these  organized  structures  etnd  their  Interactions  appear 
to  play  a  dominant  role  in  the  transport  of  heat,  mass,  and 
momentum  in  combustor  flowflelds,  the  conventional  averaging 
framework  has  not  taken  them  into  account.  Indeed,  the  concept 
of  quasi -ordered,  large-scale  structures  remains  a  challenging 
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and  hitherto  unsolved  problem  for  predictors,  since  these 
structures  are  neither  deterministic  nor  stochastic. 

3 .  PREVIOUS  WORK 

A  direct  computation  of  the  large-scale  structures  by  means 
of  the  time-dependent  Navler-Stokes  equations  has  become  an 
increasingly  viable  alternative.  The  application  of  such  an 
approach  to  the  POSF  centerbody  combustor  was  the  subject  of  the 
recently  completed  CFD  research  program  (Reference  4)  at  the 
University  of  Dayton  Research  Institute  (UDRI).  This  idea  was 
especially  desirable  in  view  of  the  photographic  evidence 
(Reference  5)  from  combusting  flow  in  the  POSF  configuration  for 
the  existence  of  nearly  ^oclsymmetric ,  toroidal  vortices 
propagating  downstream  of  the  centerbody.  The  CFD  examination  of 
the  POSF  configuration  for  nonreacting  flow  (with  only  the 
annular  air  stream  present)  involved  the  application  of  the 
computational  procedure  (Reference  6)  developed  at  AFWAL's  Flight 
Dynamics  Laboratory  (FIMM) . 

The  preliminary  results  (Reference  7)  of  the  time-dependent 
computations  and  a  computer-generated  flowfield  movie  obtained 
therefrom  appeared  to  demonstrate  the  ability  of  the  FIMM 
procedure  to  simulate  the  processes  of  formation,  growth, 
coalescence,  and  decay  of  the  wake-vortex  train.  A  comprehensive 
investigation  (Reference  8),  however,  of  the  type  as  well  as  the 
parametric  variation  of  the  boundary  conditions  at  the  inlet, 
exit,  and  solid  wall  (both  the  centerbody  and  duct  walls)  failed 
to  establish  the  previously  seen  (Reference  7)  features  of  vortex 
shedding  in  the  near  wake.  The  computations  with  the  time- 
dependent  equations  resulted  in  either  a  stationary  vortex 
[reminiscent  of  the  result  (e.g..  Reference  9)  obtained  with  the 
time-averaged  equations] ,  or  a  numerically  unstable  unsteady 
flow.  This  demonstrated  lack  of  a  numerically  stable  unsteady 
flow  from  the  time-dependent  calculations  was  of  major  concern, 
in  view  of  (a)  the  inherently  unsteady  nature  of  the  separated 


shear  layer  and  of  the  recirculating  flow  in  a  bluff-body  wake  in 
general  and  (b)  the  experimental  evidence  from  the  POSF 
configuration  in  particular. 

An  assessment  of  the  time-dependent  calculations  (Reference 
8)  emphasized  that  the  shedding-like  behavior  observed  in 
Reference  7  may  be  neither  self-sustaining  nor  physical  and  that 
a  conclusive  demonstration  of  the  successful  simulation  of  the 
dynamic  features  of  subsonic  internal  flows  was  still 
unavailable.  A  posteriori  reflection  indicates  that  such  an 
assessment  should  not  be  surprising. 

A  turbulent  flowfleld  (such  as  that  in  the  POSF  combustor)  is 

characterized  by  a  large  number  of  randomly  interacting  length 

scales  which  range  from  scales  as  large  as  a  typical  macroscopic 

dimension  (e.g.,  the  centerbody  diameter)  to  ones  as  small  as  the 

Kolmogorov  dissipative  scales.  A  computational  grid  resolving 

the  smallest  scales  would  require  grid  points  N  on  the  order  of 

Re  '  (Reference  10) ,  where  Re  is  the  representative  local 

Reynolds  number.  Estimated  values  of  N  for  realistic  flowfields 

1 3 

are  extremely  large,  on  the  order  of  10  and  higher  (Reference 
11).  Thus,  the  direct  computation  at  high  Reynolds  numbers  with 
the  full  three-dimensional,  time-dependent  equations  entails 
enormous  computing  effort.  Indeed,  with  the  present-day 
computers  and  even  those  of  the  foreseeable  future,  the 
explicitly  resolved  smallest  scale  is  very  much  larger  than  the 
dissipation  scale,  and  accurate  solutions  of  the  full  equations 
can  be  computed  only  at  very  low  Reynolds  numbers. 

Current  approach  to  direct  time-dependent  computation, 
therefore,  can  be  only  a  "large-eddy  simulation"  (LES)  at  best, 
wherein  the  full  three-dimensional  equations  are  solved  only  for 
scales  larger  than  a  selected  cutoff  length.  An  accurate  LES 
then  clearly  requires  a  "subgrld-scale"  turbulence  closure  model 
to  account  for  the  smaller  scales  below  the  LES  cutoff  length. 
While  the  selected  cutoff  between  large  and  small  scales  is 
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arbitrary  and  a  function  of  the  available  computing  power,  and 
the  fraction  of  the  spectrum  of  energy  containing  turbulent 
motions  that  can  be  explicitly  computed  will  grow  toward  unity 
as  machines  become  larger  and  faster,  this  techno logy- Imposed 
cutoff  is  unlikely  to  approach  the  order  of  the  Kolmogorov  scale 
in  high-Reynolds-number  flows. 

The  time-dependent  computations  of  the  previous  study 
(Reference  4)  are  not  truly  representative  of  LES.  A  minor 
drawback  of  the  original  FIMM  procedure  (Reference  6)  was  that 
the  calculations  were  axisymmetric,  whereas  the  large  eddies  are 
three-dimensional  and  anisotropic.  A  serious  defect  in  the  POSF 
combustor  results  obtained  in  Reference  7  and  in  subsequent 
investigations  (Reference  8)  was  the  absence  of  a  turbulence 
model.  While  the  FIMM  procedure  was  essentially  a  Reynolds- 
averaged  formulation  [with  the  closure  obtained  at  all  scales 
through  the  Baldwin-Lomax  model  (Reference  12)  of  the  mixing- 
length  type],  the  previous  POSF  results  were  obtained  from 
laminar-like  calculations,  since  no  eddy  viscosity  was  employed. 
Although  it  is  instructive  to  regard  these  calculations  (in 
References  7,8)  employing  the  molecular  viscosity  as  a  special 
subcase  of  the  Reynolds-averaged  simulation  with  zero  turbulence, 
its  relevance  to  the  high-Reynolds-number  flows  under 
consideration  is  highly  questionable.  The  argument  that  the 
numerical  procedure  does  Introduce  an  artificial  viscosity  is  not 
entirely  satisfying,  since  the  extent  to  which  the  numerical 
artifact  mimics  physical  turbulence  is  not  known. 

A  more  serious  problem  in  the  time-dependent  computations  may 
be  due  to  the  MacCormack  algorithm.  This  algorithm  was  rejected 
by  the  Stanford  University  LES  proponents,  because  of  its 
propensity  to  be  highly  energy  dissipative  (Reference  13).  Thus, 
the  asymptotic  tendency  to  a  stationary  vortex  in  the  near  wake 
observed  in  several  of  our  POSF  numerical  experiments  (Reference 
8)  may  well  have  been  caused  by  this  apparent  energy  dissipation 
of  numerical  origin.  This  conjecture  appears  to  be  supported  by 


the  analysis  of  Felerelsen,  et  al.  (Reference  13)  which  reveals 
that  when  an  initial  flowfleld  containing  all  length  scales  is 
advanced  in  time,  the  smallest  scales  (l.e.,  the  highest 
wavenumbers)  will  swiftly  disappear  through  the  highly 
dissipative  action  of  the  numerical  algorithm.  A  comprehensive 
assessment  of  time-dependent  calculations  must  address  this 
aspect . 

4.  SCOPE  OF  PRESENT  WORK 

An  assessment  of  time-dependent  calculations  for  gas  turbine 
combustor- type  flows  entails  both  unsteady-  and  steady-state 
computations.  A  comprehensive  CFD  research  of  unsteady  modeling 
for  the  class  of  ducted,  internal  flows  of  high  Reynolds  number 
and  low  Mach  number  representative  of  gas  turbine  combustor- type 
configurations  must  Involve  a  survey,  selection,  and  review  of 
potential  computational  schemes  in  the  literature,  and  an 
evaluation  of  several  time-dependent  procedures  of  relevance 
through  comparative  calculations  of  appropriate  configurations. 
The  choice  of  these  configurations  must  be  governed  by  several 
factors  such  as  geometric  simplicity,  ease  of  specification  of 
boundary  conditions  that  can  be  prescribed  or  measured,  and 
availability  of  accurate  time-averaged  and  time-resolved 
experimental  data  in  both  nonreacting  and  reacting  flows. 
Potential  candidates  which  can  satisfy  the  requirements  for  the 
test  configurations  are  likely  to  be  few.  Axisymmetric 
and  planar  mixing  layers  of  the  type  investigated  at  the 
California  Institute  of  Technology  appear  to  be  attractive 
possibilities.  The  domlnamt  interest  here,  however,  is  concerned 
with  the  POSF  research  combustor,  involving  both  centerbody  and 
noncenterbody  test  sections.  It  is  essential  that  the  evaluation 
of  the  computational  schemes  must  be  based  on  their  numerical 
accuracy  and  computational  efficiency;  their  suitability  for  use 
in  complex  combusting  geometries;  their  ability  to  be  extended  to 
include  combustion,  swirl,  and  three-dimensional  flowfields;  and 
their  accuracy  in  comparison  with  experimental  data.  Of  the 
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different  computational  procedures  available,  only  the  PIMM  code 
(Reference  6)  Involving  the  MacCormack  explicit  algorithm  has 
been  previously  examined  (Reference  8)  In  detail  for  the  POSP 
configuration.  Despite  the  major  uncertainties  in  the  MacCormack 
procedure,  alluded  to  in  Paragraph  1.3,  this  algorithm  was 
expected  to  remain  as  a  benchmark  against  which  other  com¬ 
putational  procedures  could  be  Judged. 

The  limited  scope  of  the  Scholarly  Research  Program,  however, 
meant  a  much  less  comprehensive  CPD  examination  of  the  unsteady 
flow  In  the  existing  POSP  centerbody  combustor.  Such  a  limited 
investigation,  which  is  documented  here,  essentially  involved  an 
extension  of  our  earlier  studies  (Reference  8)  by  properly 
incorporating  the  k-e  turbulence  model  of  Jones  and  Launder  (see 
References  14-16)  into  the  PIMM  (Reference  6)  calculation 
procedure  employed  previously  (References  7,8).  The  two 
additional  equations  for  k  and  £  are  solved  in  the  spirit  of 
MacCormack 's  algorithm  involving  a  forward-predictor, 
backward-corrector  sequence  at  each  time  step.  These  refined 
calculations  with  the  time-dependent  formulation  show  that  the 
flowfield  reaches  a  steady  state  asymptotically.  Purthermore, 
the  characteristics  of  the  recirculation  zone  and  the  values  of 
the  flow  variables  obtained  in  these  calculations  closely 
correspond  to  the  experimental  observations.  Present  results 
further  demonstrate  that  the  time-dependent  computations  lacking 
a  turbulence  model  for  properly  accounting  for  the  dissipation 
that  is  inherently  present  at  the  hlgh-Reynolds-number  flows  in 
configurations  such  as  the  POSP  combustor  will  invariably  lead  to 
solutions  (Reference  7)  that  are  physically  unrealistic  at  best. 

The  steady-state  calculations  were  expected  to  involve  a 
refined  version  (Reference  17)  of  the  Reynolds-averaged 
formulation.  Because  of  the  improved  physical  and  numerical 
modeling  aspects  of  Reference  17,  this  steady-state  procedure  was 
anticipated  to  yield  better  predictions  than  the  results  obtained 
earlier  (References  8,9).  Of  particular  Interest  in  this  context 
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is  the  new  POSF  combustor  configuration  under  consideration  for 
development  and  testing.  This  baseline  configuration  involving 
the  mixing  of  two  annular  air  streams  (an  outer  swirl-free  stream 
and  an  inner  swirling  stream)  and  a  central  fuel  jet  in  the  near 
wake  of  a  ducted  bluff  body  Is  expected  to  simulate  more 
realistically  the  primary  zone  of  a  gas-turbine  combustor  than 
did  the  earlier  centerbody  configuration.  Time-averaged 
computational  experiments  prior  to  and  during  the  actual  POSF 
experiments  could  serve  in  guiding  the  selection  of  optimum 
conditions  for  further  development.  The  time-averaged 
calculations  of  the  proposed  POSF  combustor  are  discussed  in 
Volume  II. 

5.  OUTLINE  OF  REPORT 

The  Incorporation  of  the  k-e  model  into  the  time-dependent 
formulation  is  presented  in  Section  II.  Section  III  reports  the 
results  of  the  time-dependent  calculations  and  shows  the 
influence  of  the  explicit  introduction  of  turbulence  in  these 
calculations.  Also  presented  in  this  section  are  some  results  of 
laminar-like  calculations  but  with  forcing  achieved  through 
time -dependent  Inflow  conditions.  The  conclusions  and 
recommendations  from  this  study  are  outlined  in  Section  IV. 


SECTION  II 

TIME-DEPENDENT  FORMULATION 

This  section  discusses  the  introduction  of  the  two-equation 
turbulence  model  for  the  calculation  of  the  nonreacting  flow  in 
the  POSF  centerbody  combustor  with  the  time-dependent  Navier- 
Stokes  equations.  The  adaptation  of  the  FIMM  computational 
procedure  (Reference  6}  is  extensively  discussed  in  Reference  8 
and  the  reader  is  referred  to  that  discussion  for  details. 

1 .  GOVERNING  EQUATIONS 

As  in  the  previous  studies  (References  7,8),  the  unsteady 
flowfleld  in  the  centerbody  configuration  is  examined  here  only 
for  the  nonreacting  situation  arising  from  the  presence  of 
annular  air  stream  alone.  Thus,  both  the  previous  and  present 
CFD  research  does  not  address  the  observed  (Reference  5)  dynamic 
behavior  in  the  centerbody  near  wake  which  involved  combustion  of 
the  central  fuel  jet  with  annular  air  stream. 

The  set  of  governing  equations  comprising  the  time-dependent, 
axlsymmetric ,  compressible  Navier-Stokes  equations  of  a  perfect 
gas  and  describing  the  conservation  of  mass,  momentum,  and  energy 
may  be  written  in  conservative  form  in  physical  space  as  follows: 

at  at  r  ar 

Here,  the  column  vectors  E,  F,  G,  and  H  respectively  represent  the 
fluxes  corresponding  to  the  time-dependent  terms,  axial  and 
radial  convective  terms,  and  the  source  terms.  These  are 
explicitly  presented  as: 
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In  Equations  (2)  through  (4)  u  Is  ths  effsctlve  viscosity 
given  by: 

M  -  M|+  (5) 

where  Is  the  solecular  viscosity  (which  is  specified  according 
to  Sutherland's  law),  and  is  the  turbulent  eddy  viscosity. 

The  latter  is  obtained  froa: 

Mt  ■  c^pk*/€,  (6) 

where  c^  is  usually  taken  to  be  a  constant  equal  to  0*09.  A  is 
the  effective  second  viscosity  coefficient  and  is  the 
corresponding  turbulent  eddy  coefflcent.  It  is  assumed  that 
A  -  -  (2/3)|i  and  "  "  (2/3)m^.  k  is  the  effective  coefficient 
of  thermal  conductivity  given  by: 
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where  the  laminar  and  turbulent  Prandtl  numbers  Pr^  and  Pr^  are 
taken  as  0*72  and  1  respectively,  and  Is  the  constant-pressure 
heat  capacity.  The  remaining  parameters  and  constants  In  the 
turbulence  model  are  =  1,  =  1*3,  =  1»44  and  Cg  *  1*92. 

It  must  be  noted  that  except  for  the  addition  of  the 
differential  equations  for  k  and  e  and  the  corresponding 
expressions,  parameters,  and  constants  for  the  k-e  turbulence 
model,  the  set  of  governing  equations  Is  the  same  as  the  one  used 
In  the  l€uiilnar-llke  calculations  of  References  7  and  8.  It  Is 
assumed  here  that  the  turbulent  gas  motion  Is  described  by  the 
two-parameter  (k-e)  model  of  Jones  and  Launder  (see  References 
14-16),  but  with  the  nonstationary  terms  taken  Into  account. 
Additional  discussion  of  the  k-e  model  and  the  refinements 
therein  to  account  for  the  streamline  curvature  and  the 
preferential  Influence  of  normal  stresses  In  the  dissipation 
equation  (these  effects  were  not  considered  for  the  present 
study)  Is  available  In  References  8  and  9. 

Finally,  the  pressure  Is  related  to  temperature  by  the 
equation  of  state 


p  =  pRT.  (8) 

The  governing  differential  equations  are  advanced  In  time  by 
MacCormack's  explicit  and  unsplit  algorithm.  The  FIMM  procedure 
of  Shang  (Reference  6)  Is  further  modified  to  solve  for  the 
additional  equations  for  k  and  e.  These  two  equations  are  solved 
In  the  same  spirit  as  MacCormack's  algorithm  Involving  a  forward- 
predictor,  backward-corrector  sequence  at  each  time  step. 

2.  BOUNDARY  AND  INITIAL  CONDITIONS 

The  boundary  conditions  for  the  inflow  and  outflow  boundaries 
and  at  the  solid  wall  are  the  following: 
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Inflow  Boundary  Conditions 


a . 

p,  pU,  pV:  Specified 


k  =  0*03U®  (9) 

and 

e  =  c  k’  •  V*  , 

where  *  =  0*03  {R(j-Rc)  .  Rd  being  the  respective  radii  of 

the  duct  and  centerbody. 

b .  Outflow  Boundary  Conditions 

The  boundary  conditions  for  the  exit  boundary  are  as 

follows : 


P  =  P, 


and 


3»ti 

3x 


0, 


(10) 


where  <^  =  U,  V,  T,  k  and  e. 

c .  Solid-Wall  Boundary  Conditions 

No-slip  conditions  (U  =  0  and  V  =  0)  are  used  along  the 
duct  and  centerbody  walls  to  specify  the  velocity  components. 

The  wall  temperature  is  a  specified  constant.  The  wall  pressure 
is  calculated  from  the  boundary- layer  assumption  that  its  normal 
derivative  vanishes  at  the  wall.  The  density  is  calculated  from 
the  equation  of  state.  The  two  scalars  of  the  turbulence  model 
are  given  by  k  =  0  and  9e/an  =  0  (where  n  is  the  direction  normal 
to  the  wall ) . 
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The  foregoing  boundary  conditions  consist  of  the  appropriate 
conditions  for  the  lamlnar-llke,  time-dependent  calculations  of 
Reference  8  and  those  of  k  and  e  equations  for  the  Reynolds- 
averaged,  time- Independent  calculations  of  References  8  and  9. 

3 .  COMPUTATIONAL  DETAILS 

The  computational  grid  employed  for  these  calculations  is  the 
baseline  domain  consisting  of  60  axial  nodes  and  46  radial  nodes. 

Figure  1  shows  this  finite-difference  grid  which  was  used  in  the 

laminar- like  calculations  of  Reference  7  and  in  some  cases  of 

similar  calculations  of  Reference  8.  Additional  details  of  this  i 

grid  are  found  in  the  latter.  Note  that  exponential  grid 

stretching  is  used  to  adequately  resolve  the  flowfield  in  the 

anticipated  regions  of  large  gradients.  Also,  for  the  present 

low  subsonic  flows,  an  artificial  viscosity  (other  than  that 

Inherent  in  the  MacCormack  scheme  itself)  is  not  needed  due  to 

the  lack  therein  of  the  very  strong  flowfield  gradients  typically 

present  in  supersonic  flows. 

Although  the  emphasis  of  the  present  study  was  the  explicit 
Introduction  of  a  realistic  turbulence  model  into  the  time- 
dependent  formulation  and  the  resulting  flowfield  behavior  of  the 
POSF  configuration,  additional  time-dependent  calculations  were 
made  to  examine  the  effect  of  external  forcing.  Both 
laminar-like  and  the  turbulent  situations  were  considered. 

a.  Laminar-Like  Cases  with  Time-Dependent  Inflow 

These  calculations  used  the  extended  (80  x  46)  finite- 
difference  grid  of  Reference  8.  Three  different  cases  were 
studied,  with  the  combustor  inlet  mass  flow  being  2  kg/s  for  all 
of  them.  A  sinusoidal  time-dependent  perturbation  was  Imposed  on 
one  of  the  inflow  boundary  conditions,  with  the  frequency  being 
128  Hz  corresponding  to  the  fiindamental  frequency  in  the 
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Difference  Mesh  (60  x  46). 


mass- flow  variations  observed  for  the  earlier  computations  (see 
p.  153  in  Reference  6).  The  unperturbed  parameters  were  akin  to 
those  of  Case  7a  In  Reference  8. 


(1) 

as  follows: 


(2) 

following: 


Case  1 

The  boundary  conditions  employed  in  this  case  are 


(a)  Inflow  Conditions 

p  =  1*1965  kg/m* 

pU  =  p48*84  (1  +  0*1  sin  ft) 

pV  =  0 


(b)  Outflow  Conditions 

p  =  9,821*39  kg^/m*,  r  <  R^ 

I®  =  0,  ^  *^d 

3U  8V  aT 

a  —  B  s  0 

dx  ax  ax 

Case  2 

The  inflow  and  outflow  conditions  were  the 


(a)  Inflow  Conditions 


p  =  1*1965  kg/m* 


pU  =  p48*84[l  +  0*1  (sin  ft  +  0*4  sin  3ft)] 


pV 

ax 


0 

0 


(b)  Outflow  Conditions 

p  *»  9,821*39  kg^/m* 

ly  =  iY  .  ar  ^ 
ax  “ax  dx  “ 
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(3 )  Case  3 


The  following  Inflow  and  outflow  conditions  were 

employed : 

(a)  Inflow  Conditions 

p  *  1«1965  kg/m* 
pU  =  p48«84 
pV  =  0 

p  =  p(2,  J)  (1  +  0‘1  sin  ft) 

(b)  Outflow  Conditions 


p  =  9,821*39  kg^/m* 

ly  =  iV  a  H  =  o 

b.  Turbulent  Case  with  Time-Dependent  Inflow 

Only  one  case  of  time -dependent  calculations  with  the  k-e 
model  was  considered  with  the  forced  inflow  condition.  The 
perturbation  considered  here  is  the  same  as  that  in  Case  2  of  the 
laminar- like  calculations  [see  Paragraph  II.3.a(2a)]. 
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SECTION  III 

RESULTS  AND  DISCUSSION 


This  section  presents  the  results  of  the  time-dependent 
calculations  with  the  k-e  turbulence  model  and  those  with  the 
time-dependent  inflow  conditions  for  both  laminar- like  and 
turbulent  cases.  These  numerical  calculations  were  performed  on 
the  NASA  Lewis  Research  Center  CRAY-IS  computer. 

1 .  INFLUENCE  OF  TURBULENCE  MODEL 

Figure  2  shows  the  velocity-vector  plots  and  vorticity- 
contour  plots  corresponding  to  the  flowfield  at  10,000,  20,000, 
30,000,  and  40,000  time  steps.  These  results  clearly  show  the 
attainment  of  a  steady-state  recirculating  flow  in  the  bluff-body 
near  wake.  It  is  interesting  to  note  that  even  by  the  first 
10,000  time  steps,  the  initial  transients  have  been  washed  off. 
Furthermore,  recall  our  earlier  observations  in  Reference  8 
(e.g..  Figures  58  and  59),  according  to  which  (a)  the  only 
unsteady  flowfield  feature  was  the  slow  axial  propagation  of  the 
recirculation  region  toward  the  exit  boundary  (at  an  approximate 
speed  of  3  ra/s) ,  and  (b)  the  appearance  of  the  reverse  flow  at 
this  boundary  triggered  mass-flow  oscillations.  With  the 
introduction  of  the  k-e  turbulence  model,  this  slow  stretching 
and  extension  of  the  recirculation  vortex  appear  to  have  been 
eliminated.  Thus,  continuation  of  the  computations  even  up  to 
40,000  time  steps  does  not  encounter  the  downstream  migration  of 
the  recirculating  region,  the  approach  of  the  reverse  flow  to  the 
exit  boundary,  and  the  consequent  incompatibility  of  the  boundary 
conditions  therein.  These  results,  therefore,  appear  to 
substantiate  our  conjecture  in  Reference  8  that  the  lack  of  an 
adequate  turbulence  model  might  have  caused  the  apparent  unsteady 
behavior  seen  in  References  7  and  8.  Furthermore,  the 
disappearance  of  the  shedding-like  behavior  (observed  earlier 
during  the  time  period  aissoclated  with  the  mass-flow  fluctuations 
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Figure  2(a).  Velocity-Vector  Plots  for  Turbulent  Plow 


Figure  2(b).  Velocity- Vector  Plots  for  Turbulent  Plow 


p  triggered  by  the  Inadequacy  of  the  outflow  boundary  conditions  to 

handle  the  appearance  of  reverse  flow  there)  resulting  from  the 
k-e  calculations  lends  credence  to  our  claim  that  the  dynamic 
behavior  of  the  POSF  combustor  predicted  In  Reference  7  is 
neither  self-sustaining  nor  physical,  but  is  strictly  a  numerical 
artifact . 

That  the  time-dependent  solution  of  the  POSF  flowfleld 
attains  a  steady  state  asymptotically  after  40,000  time  steps 
Indicates  a  failure  of  the  time-dependent  formulation  (subject  to 
time- independent  boundary  and  initial  conditions)  to  yield  a 
time-dependent  flowfleld  in  the  Interior.  This,  of  course,  is 
neither  surprising  nor  undesirable,  since  the  use  of  the 
time-dependent  formulation  to  obtain  the  steady-state  solution 
asymptotically  for  large  times,  when  the  physical  problem  does, 
in  fact,  possess  a  steady  solution,  is  well  known  in  CFD. 

Indeed,  such  an  approach  has  been  profitably  exploited  by 
Drummond  (Reference  18)  in  his  recent  calculations  of  the 
turbulent  reacting  flow  in  a  ramjet  dump  combustor.  These 
calculations  Involved  the  unsplit  MacCormack  algorithm  and  a 
mixing-length  type  algebraic  eddy  viscosity  model  (Reference  19) 
to  integrate  the  governing  equations  until  a  steady-state 
solution  is  reached.  The  calculated  results  were  seen  to 
indicate  fair  agreement  overall  with  the  experimental  results  for 
nonreacting  flowfleld.  The  observed  overprediction  of  the  rate 
of  mixing  of  air  and  hydrogen  streams  is  attributed  in  Reference 
18  to  the  use  of  the  simple  algebraic  turbulence  model  (Reference 
19)  which  is  not  considered  adequate  to  model  the  very  complex 
turbulent  flow  of  this  configuration.  The  key  point  of  Interest 
to  our  present  study  is  that  the  use  of  an  even  relatively  crude 
mixing-length  turbulence  model  in  the  time-dependent  formulation 
has  given  rise  to  a  steady-state  solution  which  is  consistent 
with  experimental  data. 

This  raises  the  Interesting  question  as  to  how  the  steady- 
state  solution  from  the  time-dependent  formulation  compares  with 
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the  time-averaged  solution  and  experimental  data  in  the  POSP 
configuration  (References  8,9).  A  conclusive  answer  to  this 
question  requires  an  examination  of  the  computed  results  of  the 
velocity  field  (corresponding  to  the  flowfield  of  Figure  2)  for 
both  the  mean  and  rms  quantities.  Such  an  examination  should 
address  the  axial  (centerline)  variation  of  the  mean  and  rms 
axial  velocity  components  (with  the  assumption  of  isotropy 
invoked  necessarily  for  the  latter) ,  as  well  as  the  radial 
distribution  at  different  axial  stations  of  the  axial  and  radial 
velocity  components  (for  both  mean  and  rms  quantities).  Such  an 
examination  was  not  within  the  scope  of  the  present  study. 

An  examination  of  the  velocity-vector  plots  and  the 
vorticity-contour  plots  of  Figure  2,  however,  offers  some 
clues.  It  is  clearly  seen  that  the  centerline  rear  stagnation 
point  in  the  velocity-vector  plot  occurs  at  a  distance  of 
approximately  one  centerbody  diameter  downstream  of  the 
centerbody.  This  result  is  consistent  with  both  experimental 
results  and  time-averaged  calculations  (References  8,9).  Another 
aspect  demonstrating  fair  comparison  relates  to  the  vortex  center 
of  the  large  stationary  vortex  behind  the  centerbody.  The  axial 
and  radial  coordinates  of  the  vortex  center  (where  both  the  axial 
and  radial  mean  velocity  components  vanish)  in  Figure  2d  exhibit 
fair  agreement  with  those  obtained  by  the  t ime- independent , 
Reynolds-averaged  calculations  (References  8  and  20) .  Note  that 
while  the  calculations  in  Reference  8  were  based  upon  an  earlier 
version  of  the  TEACH  procedure  (Reference  9),  the  results  in 
Reference  20  were  obtained  from  a  more  refined  version  (Reference 
17)  of  the  TEACH  program.  Thus,  there  is  some  evidence  that  the 
steady-state  results  obtained  here  from  the  time-dependent 
calculations  with  the  k-e  model  reasonably  conform  to  the  results 
of  time-averaged  measurements  and  calculations,  thereby  leading 
to  an  internal  consistency  between  the  time-averaged  behavior 
and  the  steady-state  behavior  attained  with  the  time-dependent 
formulation.  This  strongly  suggests  that  a  successful  CFD 
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demonstration  of  the  capability  of  the  time-dependent  formulation 
to  predict  in  nonreacting  POSF  flowfleld  the  dynamic  features 
observed  in  combusting  flow  (Reference  5)  remains  to  be  carried 
out . 

In  this  context,  it  is  of  Interest  to  point  out  that  the 
experimental  observations  (Reference  5)  of  the  flame  balls 
propagating  downstream  of  the  POSF  centerbody,  which  revealed  the 
presence  of  successive  flame  bursts  and  relative  quiescence  of  a 
random  nature,  may  arise  from  an  altogether  different  mechanism. 
This  mechanism  is  the  coupling  between  combustion  heat  release 
and  duct  acoustics.  Indeed,  because  of  the  nonlinear  coupling 
between  the  governing  equations  for  the  conservation  of  momentum 
and  energy,  confined  combustion  is  basically  oscillatory.  Such 
forced  oscillations  in  combustion  and  their  connection  with  the 
Rayleigh  criterion  which  governs  the  coupling  between  the 
mechanical  energy  in  the  pressure  wave  and  the  thermal  energy 
available  from  combustion  heat  release  are  well  known  in  ramjets, 
turbojet  afterburners,  industrial  furnaces,  and  rocket  engines. 
If,  in  fact,  the  observed  POSF  dynamic  features  (Reference  5)  are 
characteristic  of  a  forced  combustion  system,  the  rationale 
(Reference  7)  for  seeking  the  dynamic  behavior  in  nonreacting 
situations  through  a  time-dependent  formulation  becomes  suspect, 
since  the  essential  coupling  between  duct  acoustics  and 
combustion  heat  release  is  missing  in  such  a  formulation.  Here 
we  recall  the  heuristic  arguments  presented  by  Reference  21  on 
the  basis  of  time-averaged  calculations  and  perturbations  thereof 
to  hint  at  the  possibility  of  forced  oscillations  in  the  POSF 
combustor . 

Conjectures  on  the  behavior  of  the  POSF  configuration  as 
representing  a  forced  combustion  system  become  more  persuasive 
on  different  grounds  as  well.  Available  experimental  amd 
theoretical  evidence  (e.g..  References  22-24)  on  the  bluff-body 
flame  stabilizers  suggests  that  the  flowfleld  structure  during 
combustion  and  heat  release  exhibits  drastic  changes  as  compared 


to  an  isothermal  stream.  The  confined  flowfield  behind  a 
two-dimensional  bluff  body  has  been  observed  (Reference  22)  to 
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become  stationary  at  Reynolds  numbers  on  the  order  of  10  ,  due  to 
the  introduction  of  combustion,  in  contrast  to  the  flowfield  that 
is  characterized  by  large-scale  unsteady  motion  observed  in 
nonreacting  flowflelds.  Such  a  transition  from  a  time-dependent 
flow  to  a  time- Independent  flow  caused  by  combustion  has  also 
been  verified  (Reference  23)  by  a  Reynolds-averaged  formulation 
(with  the  Ic-e  model)  which  retains  explicit  time  dependence  for 
the  two-dimensional  flow  downstream  of  a  ducted  bluff  body. 
Finally,  the  recent  experiments  (Reference  24)  on  the  unconfined 
analog  of  the  POSF  centerbody  configuration  (albeit  at  small 
enough  Reynolds  numbers  for  the  flow  to  be  laminar)  have  also 
indicated  the  presence  of  vortex  shedding  from  the  bluff  body  in 
a  nearly  zero-heat-release,  isothermal  reacting  flowfield  (for  a 
TiCl^  -  H^O  reaction) ,  and  the  lack  thereof  in  -  air 

combustion  experiments. 

The  foregoing  makes  it  clear  that  what  is  not  open  to  doubt 
is  the  existence  of  the  Inherent  unsteadiness  of  a  separated 
shear  layer  downstream  of  the  trailing  blunt  end  of  the 
centerbody  and  the  initial  appearance  of  the  Kelvin-Helmholtz 
instability  waves  therefrom  (note  that  very  close  to  its  origin 
where  the  displacement  thickness  is  very  small  compared  to  the 
radius,  this  separated  shear  layer  resembles  the  plane  mixing 
layer  downstream  of  a  splitter  plate  and  is  extremely  unstable  to 
infinitesimal  disturbances).  What  is  beyond  doubt  is  that  in 
both  axisymmetrlc  and  two-dimensional  configurations,  properly 
performed  calculations  of  the  time-dependent  Navler-Stokes 
equations  with  a  turbulence  model  (irrespective  of  its  crudeness) 
do  not  capture  the  dynamic  features  but  asymptotically  lead  to 
stationary  states  for  large  times.  Therefore,  the  dynamic 
behavior  shown  in  Reference  7  is  not  correct.  This  conclusion  is 
further  strengthened  by  our  calculations  with  time-dependent 
Inflow  conditions. 
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2.  EFFECT  OF  TIME-DEPENDENT  INFLOW  CONDITIONS 

The  results  of  time-dependent  calculations  to  investigate  the 
effect  of  external  forcing  are  examined  here,  first  for  laminar- 
like  cases  (see  Paragraph  II. 3. a),  and  then  for  one  set  of 
turbulent-flow  calculations  (see  Paragraph  II. 3. b). 

a.  Laminar-Like  Calculations 

The  velocity-vector  plots  and  the  vorticlty-contour 
plots  at  different  time  steps  for  the  three  laminar  cases 
[discussed  in  Paragraphs  II.3.a(l)  through  (3)]  are  presented  in 
Figures  3  through  5  respectively.  As  noted  earlier,  these 
calculations  were  completed  for  the  extended  domain  and  thus 
permit  their  proper  comparison  with  the  corresponding  unforced 
computations  of  Reference  8.  Present  results  show  the  initial 
sheddlng-like  tendencies,  the  approach  to  a  single  recirculating 
vortex  in  the  near  wake,  and  its  slow  migration  toward  the  exit 
boundary.  The  three  different  cases  of  forcing  the  inflow 
condition  exhibit  the  same  kind  of  flowfleld  behavior  observed 
previously  (Reference  8).  There  is  no  doubt  that  with  the 
time-dependent  Inflow  perturbations,  there  will  be  mass-flow 
fluctuations  differing  between  the  inlet  and  exit  (these  are  not 
displayed  here  but  they  will  be  reminiscent  of  the  behavior 
observed  in  Figures  58  and  59  of  Reference  8),  It  is  clear, 
however,  that  the  overall  flowfield  does  not  exhibit  a  locked-on 
oscillatory  behavior.  It  is  conceivable  that  the  amplitudes  of 
the  time-dependent  perturbations  are  not  sufficiently  large  for 
this  lock-on  to  occur,  even  though  the  perturbing  frequency 
corresponds  to  the  fundamental  frequency  of  the  quarter-wave 
resonator  for  the  longitudinal  oscillations.  There  is  a  reason 
to  suspect,  however,  that  the  temporal  perturbations  at  the 
Inflow  boundary  are  damped  out  due  to  physical  and/or  numerical 
causes  (recall  the  energy-dissipative  tendency  of  the  MacCormack 


Like  Calculations  (Case  1)  with  Forcing 
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Figure  3(b).  Vortlclty-Contour  Flote  (or  Lemlnar-Llke  Oelcul.tlons  (Case  2)  with  Forcing 
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Calculations  (Case  2)  with  Forcing. 
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Figure  4(c).  Vortlclty-Contour  Plots  for  Lamlnar-Llke  Calculations  (Case  2)  with  Forcing 


Figure  4{<a)  .  Vortlclty-Contour  Plots  for  Laninar-Lllce  Calculations  (Case  2)  with  Forcing. 
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30,000 


Figure  5(b).  Vortlclty-Contour  Plots  for  Laainar-Llke  Calculations  (Case  3)  with  Forcing 


algorithm  discussed  toward  the  end  of  Paragraph  1.3).  Thus, 
while  the  present  perturbed  results  are  not  conclusive,  they  are 
further  suggestive  of  the  earlier  conclusions. 

b.  Turbulent  Flowfield  Results 

More  definitive  insights  appear  to  be  indicated  by  the 
calculations  with  the  k-e  turbulence  model.  The  inflow 
perturbation  of  this  computation  was  the  same  as  that  of  the 
second  laminar  case  [see  Paragraph  II.3.a(2a)].  As  in  the  case 
of  the  unforced  turbulent  calculations  discussed  in  Paragraph 

111.1,  the  perturbed  results  are  based  on  the  (60  x  46)  finite- 
difference  grid. 

Figure  6  shows  the  velocity-vector  plots  and  vortlclty- 
contour  plots  corresponding  to  different  times  of  the  perturbed 
turbulent-flow  computations.  While  the  vector  plot  at  60,000 
time  steps  indicates  that  the  recirculation  vortex  is  pushed 
farther  upstream  (note  the  decrease  in  the  aucial  coordinate  of 
the  vortex  center)  and  that  the  reverse  flow  is  present  along  the 
centerline  for  the  entire  extent  of  the  combustor,  at  70,000  time 
steps,  the  vortex  has  moved  downstream  (occupying  a  position 
similar  to  that  in  Figure  2b  of  the  unperturbed  flow)  and  the 
centerline  rear  stagnation  point  is  again  at  a  distance  of 
approximately  one  centerbody  diameter.  Although  detailed 
comparisons  of  the  velocity  fields  obtained  without  and  with 
forcing  would  be  helpful  in  ascertaining  the  influence  of 
time-dependent  inflow  conditions  on  the  time-dependent 
turblent-f low  calculations.  Figure  6  does  suggest  that  the 
particular  case  of  forcing  does  not  cause  an  oscillatory 
flowfield  in  the  interior.  Doubtless  the  question  of 
insufficient  perturbation  amplitude  alluded  to  in  Paragraph 

1 1 1.2.  a  remains  to  be  addressed  in  this  context  as  well. 
Nevertheless,  it  does  appear  that  a  properly  formulated 
turbulence  model  tends  to  damp  out  the  initial  transients  and 
causes  no  downstream  magnification  of  the  shear- layer 
Instability. 
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Figure  6(b).  Vortlclty-Contour  Plots  for  Turbulent  Plow  with  Forcing. 


SECTION  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

This  section  summarizes  the  main  conclusions  of  the 
assessment  of  time-dependent  calculations  for  gas  turbine 
combustor-type  flows  and  outlines  our  recommendations  for  further 
computational  fluid  dynamic  research. 

1 .  CONCLUSIONS 

The  conclusions  presented  here  stem  from  CFD  research  of 
limited  scope  addressing  the  nonreacting  flowfield  due  to  annular 
air  stream  in  the  POSF  centerbody  combustor.  No  comprehensive 
examination  of  various  time-dependent  computational  procedures 
was  undertaken.  Attention  in  the  present  study  was  limited  to 
the  FIMM  procedure  involving  the  MacCormack  explicit  algorithm 
and  focused  on  the  incorporation  of  the  k-e  turbulence  model  in 
the  time-dependent  formulation.  Thus,  this  study  represents  the 
solution  of  the  Reynolds-averaged  Navier-Stokes  equations  which 
retain  the  time-dependent  terms.  In  view  of  the  axisymmetrlc 
formulation  and  the  introduction  of  turbulence  model  for  all 
scales,  the  present  study  is  not  a  true  large-eddy  simulation 
encompassing  a  subgrid-scale  turbulence  model.  The  major 
conclusions  of  this  study  follow: 

•  Calculations  using  the  time-dependent  Navier-Stokes 
equations  with  a  properly  formulated  k-e  turbulence 
model  lead  to  a  stationary  recirculating  vortex  in  the 
near-wake  region  of  the  centerbody  for  the  POSF 
combustor . 

•  This  asymptotic  tendency  to  reach  a  steady-state  solution 
at  large  times,  which  was  suggested  by  our  earlier 
lamlnar-llke  calculations,  now  confirms  those  results  and 
clearly  demonstrates  that  the  previous  unsteady  feature 
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of  a  slow  migration  of  the  recirculating  vortex  was 
caused  by  the  lack  of  a  turbulence  model  In  the  time- 
dependent  formulation. 

•  The  asymptotically  reached  steady-state  solution  of  the 
time-dependent  problem  exhibits  a  recirculating  vortex 
which  compares  reasonably  well  with  that  obtained  from 
the  solution  of  time-averaged  equations  with  the  k-e 
turbulence  model.  This  agreement  extends  to  the 
locations  of  both  the  vortex  center  and  the  centerline 
rear  stagnation  point. 

•  Calculations  employing  time-dependent  perturbations  of 
Inflow  boundary  conditions  for  both  laminar- like  and 
turbulent  flows  do  not  show  any  oscillatory  behavior  In 
the  Interior  of  the  POSF  configuration.  However,  these 
results  with  time-dependent  bo\indary  conditions  are 
neither  comprehensive  nor  complete  enough  to  offer 
definite  conclusions  about  the  capability  to  predict 
forced  oscillations  In  the  POSF  combustor. 

•  The  attainment  of  steady-state  solutions  shows  that  the 
FIMM  procedure  does  not  capture  the  dynamic  features  of 
the  confined  recirculating  turbulent  flowfleld  of  the 
POSF  combustor.  This  was  strongly  suggested  previously 
by  the  laminar- like  computations  with  different 
combinations  of  Inflow  and  outflow  boundary  conditions  In 
two  different  finite-difference  grids  and  Is  clearly 
demonstrated  now  by  the  turbulent-flow  computations  with 
k-c  model. 

•  The  shedding- like  behavior  reported  by  the  earlier 
preliminary  calculations  Is  not  supported  by  both  the 
subsequent  laminar- like  computations  and  the  present 
turbulent-flow  calculatlona  and  must  therefore  be 
considered  suspect.  Indeed,  this  conclusion  remains 
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tenable  after  a  recent  exchange  of  views,  published 
subsequent  to  our  submission  of  the  draft  Final  Report. 
When  some  of  the  apparent  deficiencies  in  the  preliminary 
calculations  which  led  to  the  shedding- like  behavior  of 
numerical  origin  were  corrected,  the  FIMM  calculation 
procedure  failed  to  predict  the  dynamic  behavior  of  the 
POSF  flowfield.  We  believe  that  it  is  Important  to 
address  (a)  the  fact  that  the  laminar- like  computations 
with  an  extended  domain  (which  obviated  the  incompatibil¬ 
ity  of  the  exit-boundary  conditions)  led  to  a  stationary 
vortex  in  the  near  wake  and  (b)  the  corroborative  results 
reported  elsewhere  for  a  ramjet  dump  combustor  for 
which  the  time-dependent  formulation  involving  the 
MacCormack  algorithm  and  the  algebraic  mixing-length 
turbulence  model  gave  rise  to  a  steady-state  solution. 

The  key  point  of  the  present  study,  therefore,  is  to 
emphasize  that  even  relatively  crude  turbulence  models, 
when  incorporated  in  a  time-dependent  formulation  and 
properly  computed,  predict  steady-state  solutions 
asymptotically  for  large  times  and  that  the  accounting 
for  turbulence  dissipation  could  well  be  the  essential 
ingredient  that  eliminates  a  violent  self -sustained 
oscillation  unrealistically  obtained  by  a  laminar- like 
computation  of  high-Reynolds-number  flows. 

The  experimentally  observed  vortex  shedding  in  combustion 
flows  of  the  POSF  combustor  is  less  likely  to  be  an 
example  of  unsteady  flow  caused  by  the  instability  and 
its  amplification  in  the  separated  shear  layers  than  one 
due  to  the  coupling  between  combustion  heat  release  and 
duct  acoustics. 

2  See  Readers'  Forum,  AIAA  Journal  24.  April  1986,  pp.  698-701: 
M  S.  Raju,  M.  J.  Creed,  and  L.  Krlshnamurthy,  Comment  on 
"Numerical  Simulation  of  Cold  Flow  in  an  Axisymmetrlc 
Centerbody  Combustor . ” 

J.  N.  Scott  and  W.  L.  Hankey  Jr.,  Reply  by  Authors  to  M.  S. 
Raju,  M.  J.  Creed,  and  L.  Krlshnamurthy. 
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2 .  RECOMMENDATIONS 


This  limited  study  has  not  addressed  the  evaluation  of 
time-dependent  computational  schemes  from  the  viewpoint  of 
numerical  accuracy  and  computational  efficiency,  applicability  to 
complex  turbulent  flows,  extension  to  include  combustion,  swirl, 
liquid  phase  and  three-dimensional  effects,  and  comparison  with 
experimental  data.  Nevertheless,  such  a  CFD  examination  of  gas 
turbine  combustor-type  flows  is  worthwhile  and  must  be 
encouraged.  Accordingly,  we  reiterate  and  emphasize  the 

following  recommendations  offered  in  our  previous  study:  < 

•  Continue  with  the  time-averaged  formulation  (involving 
Reynolds-averaged  equations  in  nonreacting  flows  and 
Favre-averaged  equations  in  reacting  flows) ,  but 
preferably  a  three-dimensional  one  in  those  situations 
wherein  large-scale  unsteadiness  is  not  significant.  To 
properly  account  for  anisotropic  effects,  however,  direct 
solution  of  the  equations  for  the  Reynolds  stresses  must 
be  considered. 

•  Consider  the  numerical  simulation  of  those  flows,  wherein 
a  quasiperiodic  mean  flow  having  a  characteristic 
frequency  much  smaller  than  the  characteristic  frequency 
of  the  typical  large  eddies  is  encountered,  by  means  of 
the  Reynolds-averaged  (or  Favre-averaged)  formulation 
which  includes  an  explicit  time-dependence.  Initial 
consideration  of  two-dimensional  (planar)  and  axlsym- 
metric  flows  and  the  use  of  k-e  model  could  be  followed 
by  the  extension  to  three-dimensional  geometry  and 
algebraic  stress  models. 
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Explore  numerical  simulation  through  time-dependent, 
three-dimensional  Navier-Stokes  equations  as  a  large-eddy 
simulation  approach  with  subgrld-scale  turbulence 
modeling  in  nonreacting  flows  that  may  be  regarded  as 
simpler  subsets  of  the  complex  gas  turbine  combustor-type 
flows.  Viable  approaches  should  be  subsequently  extended 
to  address  reacting  flows. 
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